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Abstract 

The core collapse of a massive star in the Milky Way will produce a neutrino burst, 
intense enough to be detected by existing underground detectors. The AMANDA 
neutrino telescope located deep in the South Pole ice can detect MeV neutrinos by 
a collective rate increase in all photo-multipliers on top of dark noise. The main 
source of light comes from positrons produced in the CC-reaction of anti-electron 
neutrinos on free protons + P ^ + n. This paper describes the first supernova 
search performed on the full sets of data taken during 1997 and 1998 (215 days of 
live time) with 302 of the detector's optical modules. No candidate events resulted 
from this search. The performance of the detector is calculated, yielding a 70% 
coverage of the Galaxy with one background fake per year with 90% efficiency for 
the detector configuration under study. An upper limit at the 90% c.l. on the rate of 
stellar collapses in the Milky Way is derived, yielding 4.3 events per year. A trigger 
algorithm is presented and its performance estimated. Possible improvements of the 
detector hardware are reviewed. 



1 Introduction 

Astronomical observations have 
customarily been carried out by de- 
tecting photons, from radio- waves to 
gamma rays, with each new range of 
energy uncovered leading to fresh dis- 
coveries. Over the years, the electro- 
magnetic spectrum has been well 
covered, leaving no new frequency- 
gaps to explore. 

Neutrino astronomy is still in 
its infancy, but it is a field which 



has been growing over the past few 
decades. Besides providing a new 
probe to investigate the universe, the 
neutrino is also fundamentally differ- 
ent from the photon in that it inter- 
acts only weakly and can thus still be 
observed after passing through large 
amounts of matter. This is a mixed 
blessing, however, since the detection 
of neutrinos on Earth is made more 
difficult by the very same reason and 
requires correspondingly large de- 
tector volumes. Furthermore, neutri- 
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nos remain undeflected by magnetic 
fields. 

Until now, the only extra-terrestrial 
sources of neutrinos that have been 
observed are the Sun and the super- 
nova SN1987A. In both cases, the 
energies seen were at or below a few 
tens of MeV. The -burst emitted 
in the 1987 supernova event was de- 
tected simultaneously by the 1MB [1] 
and Kamiokande II [2,3] water detec- 
tors, a few hours ahead of its optical 
counterpart. Due to the distant lo- 
cation of the supernova in the Large 
Magellanic Cloud (~52 kpc from 
us) , only 20 neutrinos were collected 
in total. This was still enough to 
confirm that most of the energy was 
released in the form of neutrinos and 
also to vahdate the essential predic- 
tions of models describing the mech- 
anism of gravitational collapse su- 
pernovae. SN1987A proved also that 
relevant particle physics information 
can be extracted from astrophysical 
events. The data was used to set an 
upper limit on the mass of the z/g, its 
lifetime, its magnetic moment and 
the number of leptonic flavours. 

However, although the rough fea- 
tures in terms of energies, duration 
and flux were supported by the ob- 
servations, they said next to noth- 
ing about the details of the burst 
development. For this, hundreds if 
not thousands of neutrinos would 
be needed. Furthermore, both ex- 
periments were insensitive to flavors 
other than z/g, which are expected 
to carry away most of the supernova 
energy. The detection of a type II 
supernova in the Milky Way would 
provide us with a unique opportu- 
nity to study the details of the grav- 
itational collapse of a star. Several 



low-energy neutrino detectors exist 
already (LVD, Super-Kamiokande, 
SNO [4] and Baksan [5]) which will 
be able to shed light on the core col- 
lapse of a massive star, should such 
an event occur during their lifetime. 

Amanda (Antarctic Muon and 
Neutrino Detector Array [6]) is 
one of several high energy neutrino 
telescopes. Another existing detec- 
tor is NT200 in lake Baikal [7], 
whereas others are under develop- 
ment (ANTARES [8], NESTOR [9] 
and NEMO [10]). Amanda consists 
of optical modules (OMs) buried 
1500-2000 m deep in the Antarctic 
ice sheet. Each OM is made up of 
a photo-multiplier tube (PMT) en- 
closed in a pressure-resistant glass 
vessel and connected to the surface 
electronics by an electrical cable sup- 
plying power and transmitting the 
PMT signals. 

Amanda is designed for the ob- 
servation of TeV neutrino-sources, 
utilizing the large volume of trans- 
parent glacier ice available at the 
South Pole Cherenkov medium. 
In spite of the much lower neutrino 
energies of (9(10 MeV) involved in 
a burst, it has been shown [11] that 
a detector of this type could also 
be used successfully to monitor our 
Galaxy for supernova events (the de- 
scription of a method for neutrino 
telescopes using ocean water can 
be found in [12]). Since the cross- 
section for inverse /5 decay reaction 
on protons + P ^ + in the 
ice exceeds the cross-sections for the 
other neutrino flavors and targets, 
Ve events are dominating Dur- 



Note that the scattering of on "O 
is negligible at the expected energies. 
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ing the estimated ~ 10 sec duration 
of a neutrino-burst, the Cherenkov 
hght produced by the positron tracks 
will increase the counting rate of 
all the PMTs in the detector above 
their average value. This effect, when 
considered collective behavior, 

could be seen clearly even if the in- 
crease in each PMT would not be 
statistically significant. An obser- 
vation made over a time window of 
several sec could therefore provide a 
detection of a supernova, before its 
optical counterpart is observed. The 
stable and low background noise in 
Amanda (absence of '^^K and of bi- 
oluminescence in the ice) is a clear 
asset for this method. 

In this paper, the 302 OMs of the 
Amanda BIO stage completed in 
1997 are used for the analysis. In 
the next section, we summarize the 
theoretical expectations for the neu- 
trino burst which is emitted when 
a massive star undergoes gravita- 
tional core collapse. An account is 
then given of the detector hardware, 
as well as of the data collected with 
it. The analysis method is explained 
next and the results obtained from 
its application to experimental data 
are presented. A study of the detec- 
tor performance follows. Next, the 
principle of an online trigger algo- 
rithm is described. We end with a 
discussion of possible improvements 
of the detector. Early searches for su- 
pernova neutrino-burst signals with 
the Amanda detector have been 
presented in [13-15]. 



2 Theoretical preliminEiries 

When the core of a massive star 
(M > SMq) runs out of nuclear fuel, 
it collapses and ejects the outer man- 
tle in a SN explosion of type Il/Ib/Ic 
[16]. Only ~ 1% of the energy is re- 
leased in kinetic and optical form, 
whereas the remaining 99% of the 
gravitational binding energy change, 
about 3 X 10^^ ergs, is carried away 
by neutrinos [17]. 

During the process, the inner core 
reaches nuclear densities, bringing 
the collapse to an abrupt stop, and 
a shock-wave is formed. The front 
of this shock is driven outward, 
through the in-falling material, pass- 
ing through the neutrinosphere (the 
location where the material changes 
from opaque to transparent to neu- 
trinos) on its way, until a point where 
it stalls . 

During the first 10 msec, a burst 
from the neutronization process 
e" + p — >• n + Ue releases ~ 10^^ 
ergs. Although neutrinos interact 
only weakly, the densities built up 
during the collapse are so high that 
they cannot stream out. Instead, 
they are trapped and diffuse out over 
a time scale of several sec. When 
they finally reach the neutrinosphere 
they can escape, with a thermal 
spectrum which is approximately 
Fermi-Dirac [17]. The trapped neu- 
trinos are produced in pairs through 
the reaction e~ + e"*" u + u for 
all lepton fiavors. In addition, the 
UgUe pairs are also produced via 

p{n) + e~^"^^ n{p)+ Ue- They are 
produced with distinctive energies, 
because the neutrino- spheres for each 
type are located at depths with dif- 
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ferent temperatures. The and u^- 
and their anti-particles have a mean 
energy of {E) ^ 25 MeV. The i^e have 
a mean energy {E) 16 MeV and 
the corresponding value for neu- 
trinos is (E) ^11 MeV [18]. In total, 
the radiated energy is expected to 
be equally distributed over each fla- 
vor of neutrinos and anti-neutrinos 
[19]. The luminosity-profile can be 
described by a quick rise over a few 
msec and more and then falling over 
a time of C (several sec), roughly 
like an exponential with a time con- 
stant r = 3 sec [20] . During the late 
phase the neutrino luminosity possi- 
bly follows a power law with index 
~ 1 ± 0.5 [17]. The detailed form of 
the neutrino luminosity used below 
is less important than the general 
shape features and their characteris- 
tic durations [20]. 




r[kpc] 

Fig. 1. The distribution /(r) of pro- 
genitor stars in the Milky Way, located 
within a distance r from the Earth. In- 
set: the corresponding cumulative dis- 
tribution F{r) (taken from a model by 
BahcaU & Soneira [21,22]). The steep 
behavior around ~ 8 kpc is due to the 
high star density in the Galactic bulge. 

Currently, the rate of galactic 



stellar collapses is estimated to be 
one every 40 years [23], with a 
best experimental upper limit of 
1 every 5 years [5]. At least three 
underground neutrino detectors 
(Supcr-Kamiokande [24], Baksan [5], 
LVD [25]) have a sensitivity high 
enough to cover the Galaxy and be- 
yond ( the distribution of progenitor 
stars in the Milky Way is shown in 
Fig.l). 

3 Description of the detector 

The Amanda detector has been 
deployed over several years [6]. Each 
deployment phase consists of drilling 
holes in the ice and lowering strings 
of OMs connected to electrical cables 
(see Fig. 2). In this text, we will limit 
ourselves to the study of data taken 
in 1997 and 1998, with ten strings 
(302 OMs) in operation during these 
two years. 

The first four strings were de- 
ployed in 1995-1996, using coax- 
ial cable. Strings 5-10 followed in 
1996-1997, using twisted-pair ca- 
ble. Hamamatsu photo-multipliers 
(PMTs) of type R5912-02 were used 
in all OMs. The PMTs of strings 1-4 
are enclosed in spheres manufactured 
by Billings, whereas the remaining 
strings use spheres with better trans- 
parency [14], produced by Benthos 
company. The AMANDA-supernova 
detector was part of the AMANDA- 
SNMP (SuperNova-MonoPole) data- 
acquisition system (DAQ) [14], which 
operates independently from the 
AMANDA muon-trigger and DAQ. 
It consisted of custom-made CA- 
MAC modules and is read out by 
a Macintosh via a standard SCSI 
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interface. 

The idea is to continuously mea- 
sure the counting rates of aU OM 
channels separately and store that 
data for further analysis. For this, 
each channel counted the arriving 
OM-pulses synchronously into a sep- 
arate 12- bit counter during a fixed 
time interval of 500 msec which is 
synchronized by a GPS-clock. 

Correlated noise following OM- 
pulses was suppressed by a com- 
mon programmable dead-time in the 
range of 100 ns to 12.8 jj,s (which was 
set to 10 fjLs). 



AMANDA-B10 



20 MeV 
positrons 

-A 




1 meter 



Fig. 2. View of the detector and of 
the expected positron tracks produced 
in CC z7g interactions within Amanda 

The present detector configuration 

has nine additional strings, bringing 
the total number of OMs to 680, and 
the data is currently taken with an 
upgraded VME/ Linux based DAQ. 

3.1 Low- energy in Am Al<lD A 

The anti-electron neutrinos pro- 
duced in a supernova follow a Fermi- 
Dirac distribution with (Ep 



16 MeV , which leads to a peak 
value of the positron energy distri- 
bution E^+''^ - 20 MeV [11]. The 
corresponding track length is thus 
~12 cm [26-28]. The expected sig- 
nal for a SN1987A-type supernova 
has been simulated [27,28], yielding 
a predicted excess in the number of 
photo-electrons per PMT: 



A^p.e. ~ 



11 



Pice • Kflf 



Kam 



52 kpc 



kpc 



where 11 is the number of events 
detected by Kamiokande II, MKam 
is the effective mass of Kamiokande 
II and (ikpc is the distance between 
the supernova and Earth. The value 
of MKam is 2.14 ktou, times a fac- 
tor 0.8 to correct for the fact that 
Kamiokande II had an energy thresh- 
old, whereas the AMANDA supernova 
system does not (see [28]). The ice 
density, pice, is 0.924 g/cm^. In pre- 
liminary studies, the effective volume 
Vcs of one OM was found to be ap- 
proximately proportional to the ab- 
sorption length of ice and estimated 
to be 414 m^ (using an ice-model 
with 90 m absorption length) [28] . A 
subsequent detailed treatment was 
made, based on the depth-dependent 
ice properties described in [29]. This 
showed that V^s varies, depending 
on the OM location. The resulting 
values for string 1 function of 

OM number are shown in Fig.3. The 
different absorption of Benthos and 
Billings glass spheres has also been 
taken into account. Following Eq.l, 
we expect ~100 counts/OM during 
10 sec [11] for a SN1987A-like su- 
pernova at 8 kpc distance from the 
Earth. 

The summed rate over A^om opti- 
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400- . 
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100- 



5 10 15 20 

OM 

Fie. 3. Effective volume function of 
OM number in string 1. The variations 
are due to the change in ice properties 
as a function of depth. The OMs are 
numbered from their highest to lowest 
location. 

cal modules in the presence of a su- 
pernova is: 

>S'tot >S'noise ~l~ 'S'signal (2) 
i=l 

Snoise is the dark noise summed over 
all OMs and ^signal is the summed 
rate of photo-electrons resulting from 
z/g-interactions. R^"^^*^ and Rf^^^^ are 
the rates of background and signal 
counts, respectively, for each indi- 
vidual OM. If a supernova would 
occur, one should observe a sud- 
den increase in S'tot, from S'noise to 

'S'noiso ~l~ 'S'signal [1 1] • 

3.2 The noise background 

Thanks to the low temperature in 

the ice, 35°C, the PMT noise is 

strongly reduced. Furthermore, the 
ice is a low-noise medium, free from 



natural radioactivity and biolumi- 
nescence. However, in order to detect 
a small excess due to the signal, it is 
not the level of the dark-noise S'noise, 
as much as its fluctuation that has 
to be kept low. If the dark noise from 
the OMs would be purely Poisso- 
nian, the fluctuation of S'noise during 
At = 10 sec would be: 



^5„oisc = V ■ ^noisc ■ A^OM (3) 

where -Rnoise is the typical OM noise 
rate in Hz. 

In reality, the situation has been 
found to be more complex and the 
Poisson assumption had to be mod- 
ified for various reasons. Measure- 
ments in the laboratory have re- 
vealed two types of after-pulses due 
to ionized rest gases, one delayed 
by ~1.5 fj,s and the other by ~6 fis. 
This situation has been improved, by 
implementing a 10 /xs artificial dead- 
time to suppress after-pulsing [14]. 

At that stage the OM dark-noise 
follows a Gaussian distribution with 
average rates of 300 Hz for OMs in 
strings 1-4 and 1100 Hz for OMs in 
strings 5-10. The difference between 
the two groups is due to the higher 
potassium content of the spheres 
used in the second set. One of the 
processes by which /5-decaying ^^K 
increases the counting rates is by 
scintillation!^, where the primary 
electrons produce secondary pho- 
tons. This non-Poissonian process 
has been studied and is described 
in [30]. Estimates of the potassium 
content in the different spheres used 



^ Cherenkov Hght would have a signa- 
ture of high ampHtude pulses, which are 
not observed. 
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were made in the laboratory, yield- 
ing 2% for Benthos and ~ 0.1% for 
Billings spheres. The potassium con- 
tent of the PMT glass is only 0.013%; 
moreover, the PMT envelope is much 
thinner than the glass vessel. 



3.75 r 




1.75 r 



2000 4000 6000 8000 10000 12000 14000 
#counts in 1 s 

Fig. 4. F as a function of the noise, 
showing the strong correlation with the 
average of the measured rates. Note 
that the ordinate starts aX F = 1.5. 

It was observed that the vari- 
ance of the OM dark-noise is larger 
than expected. A Fano factor can be 
used to describe this deviation from 
a Poissonian behavior (see [30]), 
F = a\N)/{N), where (N) is the 
average number of counts during a 10 
sec intervals. For a Poisson process, 
F = 1. We get F ^ 2.2 for strings 1- 
4 and F fti 3 for strings 5-10. This is 
shown in Fig.4, where the behavior 
of -F as a function of the dark noise 
is fitted with a line. 

Another part of the dark noise is 
due to Cherenkov light from atmo- 
spheric muons passing nearby. Even 
though the muon flux is considerably 
attenuated by the ice-cover, this con- 
tribution is not completely negligible. 
We have estimated a noise rate of 25 
Hz due to muon light at the topmost 



OM and 12 Hz at the bottom. 

4 Data 

The data were rebinned before 
analysis, using a time window of 10 
sec for the rate calculation instead of 
the original 0.5 sec used in the DAQ. 
Henceforth, each 10 sec interval will 
be referred to as an event and all 
rates will be expressed in [Hz]. Each 
run has a typical length 14.6 hour. 

The quahty of the data depends 
on two main factors: the stability of 
each OM over long-term periods and 
external disturbances of the DAQ 
which can render a run useless. Data- 
quality varied considerably over the 
year for several OMs. In order to 
reach a state where a fixed subset of 
OMs were stable in a fixed collection 
of runs, we developed an iterative 
selection algorithm. As an input, we 
produced a binary table describing 
the quality of each OM in each run, 
a 1 representing a ^good OM in the 
corresponding data-run and a rep- 
resenting a low quality OM in that 
run. An OM is meant to be 'good in 
a run if it exhibits a Gaussian be- 
haviour with a noise rate consistent 
with its average over a year. 

The resulting matrix was then pro- 
cessed iteratively according to the 
following steps: 

1 . The sum of each row and each col- 
umn is calculated. 

2. The OM or run corresponding to 
the smallest sum found in step (1) is 
removed from the table. Here, the ba- 
sic assumption is that OMs and runs 
are of equal value for the analysis. 

3. The iteration is stopped if there 
are no O's left in the table, i.e. when 
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S,,,, [kHz] 

Fig. 5. Top: summed counting rate Stot 
function of time for the selected 
OMs. The fuh hne represents the mov- 
ing average. Bottom: the correspond- 
ing S'tot-distribution (dots) fitted with 
a Gaussian (curve). 

each OM left is stable in each run 
left. Otherwise, the procedure is re- 
peated by going back to step (1). 
Table 1 summarizes the effect that 
the cleaning has on the live-time and 
number of selected OMs. The OM 
set is almost the same in both 1997 
and 1998 and the live-time difference 
is mainly due to calibration activity 
involving laser light in 1997, which 
caused more runs to be rejected by 
the algorithm that year. Of the 224 
OMs selected in 1997 and 235 in 
1998, 52 come from the 80 OMs in 
strings 1-4 and the remainder come 
from strings 5-10. 



1997 


Runs 


Days 


OMs 


Available 


318 


185 


302 


After cleaning 


179 


103 


224 


1998 


Runs 


Days 


OMs 


Available 


299 


164 


302 


After cleaning 


205 


112 


235 



Table 1 

Summary of the effect of data-cleaning, 
both for 1997 and 1998. 



The variation of S'tot is shown for 
one 'good' run in Fig. 5 (top). Possible 
deviations from a Gaussian distribu- 
tion which are due to slow drifts of 
the noise rates over time (see Fig.5, 
bottom) can be corrected for, as de- 
scribed below. 

5 Analysis 

Though long enough to cover the 
supernova burst when synchronized 
with it, an arbitrarily located time 
window of 10 sec will (on average) 
not match the onset of a supernova 
signal (see Fig.6). If we consider the 




Time 



Fig. 6. The location of the time win- 
dow of observation and the SN burst 
are uncorrelated, resulting in a reduced 
sensitivity. 

canonical supernova neutrino-burst 
time profile mentioned in Sec.l, the 
loss in signal due to the time-window 
location is ~ 20%. 

The OM-stability is a crucial pa- 
rameter for our purposes and it can 
be influenced by several factors. 
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Some of them, such as high-voltage 
variations, arc of known origin and 
their effect can be both understood 
and estimated. Long-time trends 
are also present. E.g., due to PMT 
aging, the dark noise decreases by 
~ 2% per year. These variations are 
beyond our control and can only be 
monitored. Other disturbances are 
more irregular, such as human inter- 
ventions on the front end electronics, 
or periods of severe weather condi- 
tions at the South Pole which have 
been observed to interfere with the 
data-taking. 

In order to remove monotonic 
trends and periodical fluctuations on 
a scale of several hours and longer, 
we subtract a moving average {MA) 
from the S^ot values. The residual 
deviation RES, defined by the dif- 
ference between S'tot and MA for 
each event (i), should then yield a 
time series stationary on scales of 
the order of an hour and longer [31]: 



RES{S,,,{^) = 
Stot{i) 



MA{Stotit)) (4) 



with the average taken over a time 
window of length (n^indow + 1) bins: 



M^(5tot(^)) 



■X 



('T'window ~l~ 1) 

^window /2 

E ^tot(^+J)(5) 

window /2 

The MA-subtraction acts like a high- 
pass filter and was checked not to 
affect our SN-burst signal detection 
ability. 

For this off-line analysis, the 
sum is taken symmetrically over 
(^window + 1) time-bins around the 
studied event and will therefore be 



referred to as MAgy^a- We found that 
1000 sec was an appropriate interval, 
correcting for most of the variations 
observed in the data, as shown in 
Fig. 7. Fig.8 shows the residuals of 




10 12 14 
time [hours] 



Fig. 7. The same time series as in Fig. 5, 
after moving-average subtraction and 
shifted by {Stot)- 

the summed noise, histogrammed for 
the 215 days of live-time for 1997 
and 1998 combined. The peak of 
the distribution is very well fit by 
a Gaussian distribution. Assuming 
that the dark noise peak is Gaussian 
and that the data from the different 
OMs are uncorrelated, we sould ex- 
pect one entry per year above a value 
of 1235 Hz and 1 entry per 215 days 
for RES > 1200 Hz. 

In order to get rid of more back- 
ground events with large RES, we 
need to develop new tools using the 
more detailed data information at 
our disposal. We can distinguish be- 
tween three classes of events in the 
experimental data: 

(1) pure dark noise background, 

(2) supernova signal -|- pure dark 
noise (i.e. our signal events), 

(3) events where instrumental noise 
is present. 

We will start by assuming that 
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5 10 15 

RES [kHz] 



Fig. 8. Distribution of residual summed 
noise RES for 215 days live-time of 
1997 and 1998 data. A Gaussian fit of 
the peak of the distribution yields a 
standard deviation a = 265 Hz. 

the dark noise, as well as the ex- 
cess counts caused by the signal, 
are Gaussian and that all OMs are 
equivalent. If a supernova would be 
responsible for an event with high 
RES, we would expect the noise 
rates of all OMs to increase coher- 
ently. The first and second classes 
can easily be separated by cutting on 
the value of the summed noise resid- 
ual RES. The third class is more 
problematic, since those events can 
yield a large RES-vaAuc due to an 
increase in noise rates caused by dis- 
turbances occurring outside the de- 
tector itself. In that case, it is reason- 
able to suspect that the parts of the 
front-end detector which are affected 
are located in the control room at the 
surface of the ice and that different 
components of the electronics pick 
up noise in different ways. If such 
events are caused by some unstable 
channels, rather than by an overall 
increase of the noise in the detector. 



this should be a clear signature for a 
background event. 

With the aim to remove these fake 
supernova candidates, we consider 
the likelihood of each event. Fur- 
thermore, since the rates are nearly 
Gaussian, a function can be used 
instead of a likelihood function: 

X=^E(^^^)^ (6) 

where Ui is the measured noise of each 
of the A^oM optical modules and /li 
is its expectation value -the moving 
average of the i:th OM will be used 
as an approximation. 

The difference (n, — /Xj) is the de- 
viation of the noise from its moving 
average. A/i represents the expected 
number of counts per OM per sec, 
which are caused by the signal, and 
comes in addition to the mean dark 
noise. Typically 100 counts in 10 sec, 
or 10 Hz per OM would be caused 
by a SN1987A-type supernova at the 
center of the Milky Way. Finally, ai is 
the standard deviation for OM num- 
ber (i). The corresponding variance 
is (7? = (ToMi + A/// At, where the 
term A/i/At is added quadratically 
to the spread of the OM dark noise 
o'oMi to account for fluctuations in 
the signal, keeping in mind that 
A/x is a rate measured over At — 10 
second but given in [Hz]. For each 
analyzed run, a Gaussian flt is made 
of the noise rate of OM number i 
after the moving- average correction 
and the resulting standard deviation 
(7oMi is used. 

Note that the same function as 
in Eq.6 would be set up in the hypo- 
thetical situation where one would 
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measure a standard supernova-signal 
Nqm times using only one OM. In 
the case of classes (1) and (2) of 
events with homogeneous noise, the 
X^-value will be small. By contrast, 
that value is expected to be large for 
any Afi if the noise recorded by the 
OMs is not uniform throughout the 
detector. By cutting on x^, one can 
thus expect to reject events of class 
(3) above. 

Near the background, i.e. when 
A/x/At << (ToMi) we can solve for 
A/x, by minimizing Eq. 6 w.r.t. A//: 



^NoM 



rii 



(7) 



so that we get an actual measure- 
ment of the supernova-induced rate 
excess per OM, A/x, predicted in 
Eq. 1. With the assumption that 
all OMs are equivalent, the variable 
RES is the product Nqm x A/x. The 
standard deviation of A// is: 



„noise 



(8) 



This expression for a"^^^^ yields the 
standard deviation of the back- 
ground distribution, i.e. when A/x = 
(see Fig. 10). The characteristics of 
the individual OMs are taken into 
account in the computation of A/x 
with the different values of CTj. 

We can also introduce a module- 
dependent sensitivity in Eq.6 to 

treat properly the variation of the 
absorption-length of the ICG clS cl func- 
tion of depth [29], modifying Eqs.6 
and 8: 



2 ^^Vn,-//,-e,-A/x^2 

X = 

1=1 



„noise 



\ 



ENqm 



(9) 



(10) 



Here ti is the relative sensitivity of 
OMj, calculated as Vgg/V^r where 
V^Q = 414 m^ is the reference volume 
mentioned is section 3.1. 
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Fig. 9. Background distribution and 
expected SN signal of A// for super- 
novae at different distances, in kpc. The 
smearing of the signal due to the fixed 
10 sec observation window is taken into 
account (see Fig.6). Notice the logarith- 
mic scale in A/x. 

Fig. 9 shows the distribution of the 
A/x signal expectation for supernovae 
of type SN1987-A located at differ- 
ent distances, including the effect of 
the fixed 10 sec time-window loca- 
tion. The expected distribution in the 
absence of a supernova is also shown, 
for comparison. 

The resulting A/x distribution for 
the combined 1997 and 1998 data- 
sets can be seen in Fig. 10, where a to- 
tal live-time of 215 days has been ac- 
cumulated. A cut on x^/n.d.f. < 1.3 
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has been applied (leaving 99.9% of 
events), removing fake events which 
could still be seen in Fig. 8. With- 
out the x^/n.d.f. cut, but cutting 
on A// > 4.0 Hz, we have 1.7 back- 
ground events per week. The spread 
of the background a^^^^'^ = 0.8 Hz is 
the same for both years. This reflects 
the fact that the eleven additional 
OMs used in the 1998 set have high 
noise levels, and their spreads con- 
tribute little to reducing the total 
spread in Eq.lO. 
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Fig. 10. Distribution of Afi after appli- 
cation of cuts, for 215 days of live-time. 
The vertical line is set at the level where 
we get a rate of one background event 
per year. 90% of supernova neutrino 
bursts located at 9.8 kpc distance would 
be seen above that cut. 

Considering the background distri- 
bution of Fig. 10, a SN signal detected 
with 90% efficiency at a distance of 
8, 4 and 2 kpc would correspond to 
an excess of 7, 29 and 116 standard 
deviations (see Fig.9). The hue on 
Fig. 10 indicates the A//-level (= 4 
Hz) above which we would expect one 
fake event per year. Since the distri- 
bution in the figure is well fitted with 
a Gaussian, this number is found by 



fitting the data and integrating from 
that point and up to infinity. A cut 
on A/i > 4 Hz corresponds to a 90% 
efficiency for a SN 1987A-like super- 
nova located at a distance of 9.8 kpc. 

6 Supernova sensitivity 

The results of the analysis, along 
with numerical simulations, can be 
used to further assess the perfor- 
mance of Amanda as a supernova- 
burst detector. 

Since the fitted noise-rate excess 
per OM is centered at zero in the 
absence of any supernova and has a 
known Gaussian spread, the number 
of dark noise background events is 
easily calculated for a given cut on 

The signal that we expect to see 
is characterized by a frequency and 
a magnitude, both of which must be 
estimated to decide how to compute 
our signal-to-noisc. The more elab- 
orate estimates of the frequency of 
gravitational collapses are made us- 
ing various methods and combining 
them [23] . We will limit ourselves to 
the most conservative guess of 1 su- 
pernova/century given in [32], since 
the performance of the detector does 
not depend strongly on that number, 
as will be shown. 

The size of the signal depends on 
the distribution of progenitor stars in 
the Galaxy [22] and on the supernova 
neutrino-luminosity. The fraction of 
stars which could undergo a gravita- 
tional collapse is shown in Fig.l. We 
decided to use the measured luminos- 
ity of SN 1987A as an estimate of the 
signal-strength, since it is so far the 
only observation of an actual event. 



13 



10" 




10 1 10 10^ 



A^[Hz] 

Fig. 11. Expected A// distributions at 
(full curve) and within (dotted curve) 
a distance of 10 kpc, shown in log-log 
scale. The dominant contribution to the 
signal is given by the closest stars, which 
appear in the tail of the distributions at 
high values. The vertical lines mark 
the 90% detection efficiency lower limit. 

Simulation studies were made in 
order to estimate the signal, for which 
we used the time-profile described in 
section 2. The effect of the 10 sec 
time- window inefficiency (see Sec. 4) , 
as well as the Gaussian smearing due 
to dark-noise fluctuations, was then 
applied to a signal of strength given 
by Eq. 1. As an example, the resulting 
A/x-distribution for a supernova oc- 
curring at a distance of 10 kpc can be 
seen in Fig. 11 (full curve). If we fold 
in the probability distribution shown 
in Fig.l, we can also obtain the signal 
distribution for supernovae occurring 
within, rather then at a radius of 10 
kpc from the Earth (see Fig. 11, dot- 
ted curve). 

With these definitions of back- 
ground and signal {i.e. everything 
within a sphere of a given radius) , it is 
now possible to define a performance 
parameter which we call signal-to- 



noise S/N at a given distance. For 
this purpose, the A/i-distribution of 
progenitor stars within that distance 
is computed as well as the lower cut 
which keeps 90% of stars above it. 
The signal S above the A/x-cut is 
calculated, taking all efficiencies into 
account, as well as the estimated 
rate. Finally, the noise N from the 
dark noise background for the same 
period of time and same cut on A/x is 
readily computed from the Gaussian 
assumptions. Fig. 12 shows the re- 
sulting S/ N parameter as a function 
of distance from the Earth. At a char- 
acteristic radius, a drastic drop in 
S/N from very large values down to 
essentially zero occurs, which char- 
acterizes the detector performance. 
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Fig. 12. The S/N-ratio as a function 

of distance for strings 1-10 {S/N > 1 
within 11 kpc, i.e. 60% of the Galaxy). 
S/N = 1 corresponds to a cut at 
A/x = 5.5 Hz. Inset: the same curve seen 
in a logarithmic scale. 

We will therefore take the point at 
which S/N — 1 as Si natural measure 
of the performance of different de- 
tector configurations. Note that the 
distance that corresponds to that 
point does not depend strongly on 
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the chosen S/N value, due to the 
sharpness of the drop. For instance, a 
larger estimate for the gravitational 
collapse rate would not change that 
characteristic distance significantly. 
Indeed, using the different estimates 
found in the literature, of 1 SN ev- 
ery 100 [32], 47 [23] and 11 [33] 
years, the S/N — 1 point would be 
located at 11 (~ 60% Galaxy cover- 
age), 11.3(63%) and 11.8(66%) kpc, 
respectively. The small differences 
are due to the steepness of the curve 
(see Fig. 12). Amanda BIO is thus 
able to detect supernovae within a 
distance of 11 kpc with a S/N > 1, 
i.e. to cover ~ 60% of the stars in 
the Milky Way with one background 
fake per century. The coverage for 
detector configurations with other 
"^aT^" values is shown in Fig. 13. If one 
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Fig. 13. Percentage of the Galaxy cov- 
erage versus cr'^^'^. Note that no specific 
detector setup is assumed in this figure. 

lowers the requirement to accepting 
one background event per year and 
a supernova detection probability of 
90%, the Galactic coverage is 70%. 
Even though the additional OMs in 
strings 5-10 have worse noise char- 
acteristics, they provide a significant 



improvement. Observe however, that 
this is mainly due to the very steep 
behavior of the distribution in Fig.l 
in the range of 7-14 kpc around the 
bulge, where the star density is the 
highest. 

The non-observation of any super- 
nova signal can be used to put an up- 
per limit on the rate of gravitational 
stellar collapses: 



where rid. is the upper bound at the 
90% c.l. on the expected number of 
events, rj is the Galactic coverage 
(70%), e is the detection efficiency 
(90%) and T is the live time (215 
days). Using the procedure described 
in [34] (with no candidates observed 
and one estimated background fake 
per year), we obtain an upper limit 
on the rate of gravitational stellar 
collapses at the 90% c.l., of 4.3 events 
per year. 

7 Amanda Supernova Trigger 
Algorithm (Asta) 

The goal of the trigger is to detect 
a supernova event in real time, in 
order to send out a prompt alarm to 
the astronomical community, ahead 
of its optical observation. A future 
implementation of AsTA at the 
South Pole will make it possible to 
join SNEWS [35], a supernova early 
warning network of neutrino detec- 
tors. A first version of ASTA based 
on 11 days of data was presented in 
[36]. 

The online algorithm must be able 
to suppress background events which 
could fake a SN neutrino burst, es- 



15 



pecially those due to non- Gaussian 
variations of the noise rates. We must 
use a moving average which uses ex- 
clusively events occurring to the left, 
i.e. before the considered time-bin. 
We win refer to it as MAidt, as op- 
posed to MAsyra, with the difference 
that the sum in Eq.5 is now taken 
from {i — n^indow) to i. In order to 
accommodate for faster variations of 
the data, we define an additional vari- 
able, the derivative of MAiett for each 
event (i): 

DER{i) = (12) 
[M^left(i) -MAieft(i-l)]/Ai 

where At — 100 sec. It has been 
tested with numerical simulations 
that during normal data-taking con- 
ditions, the occurrence of a SN signal 
should not increase DER{i) signifi- 
cantly. 

Using the same OM selection as 
in the offline-analysis was found to 
be unsuitable for ASTA, rejecting a 
large percentage of the events after 
applying a cut on x^/n.d.f. < 1.3. 
Instead, we removed 117 known un- 
stable OMs and then applied the 
cleaning algorithm described in sec- 
tion 4 on the 1998 data sample. This 
resulted in a set of 170 OMs that 
were then used to test the cuts on 
both the 1997 and 1998 data. Fast 
time variations of the data (S'tot) 
are not always corrected fully by 
moving-average subtraction. Taking 
this into account we decided to it- 
erate, i.e. apply the moving-average 
correction a second time, on RES: 

RES,te,(RES{{}) = 

RES{{) - MAi^ft{RES{{)) (13) 

A new A/i can be calculated, as 
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Fig. 14. Tail of the background distri- 
bution after different cut levels: (dotted 
line) 'no cuts', (dashed line) + DER 
cut', (solid line) 'all cuts'. The vertical 
line indicates the A^u-cut that leaves one 
statistical fake per year. The live time 
is 296 days. 

well as new variables xf^^j, and 
DERiter- The final cuts applied 
are: A// > 4.1 Hz, \DER\ < Is'^, 
xVn.d.f. < 1.3, \DERitcr\ < Is'^ 
and xLr/'^-'^-f- < 1-3) leaving five 
triggers only, i.e. one every two 
months (see Fig. 14). Note that these 
requirements leave the efficiency 
to trigger on a SN during stable 
data-taking conditions unaffected. 
In order to estimate the perfor- 
mance of the cuts as previously, we 
use (7^7^ = 0.85 Hz and the cut 
A/j, > 4.1 Hz, resulting in a 65% 
coverage of the Galaxy with 90% 
detection efficiency. On top of the 
five fakes found experimentally, we 
expect statistically one background 
event per year. We note that all fakes 
have low significance and that al- 
ready a cut on A/x > 5 Hz can get 
rid of them, which would lower the 
Galactic coverage to 55%. 
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8 Discussion 

The good overall performance of 
the detector is mainly due to the 
low noise of the OMs in strings 1-4. 
However, even if all OMs would be 
of the same type as those, 800 of 
them would still be required to cover 
100% of the Milky Way Prom Eq.lO 
and Fig. 13, it is clear that a'^^ is 
a crucial parameter for the detector 
performance and must be kept as 
small as possible. 

Assuming a uniform detector with 
only one type of OMs, Eq.lO be- 
comes: 



(14) 



OM 



One sees from Eq.l4 that increasing 
A^OM is a less efficient way of reduc- 
ing cr^^'^ than increasing the collec- 
tion efficiency e or reducing aoM- the 
spread in the dark-noise rate of each 
OM. Several technical improvements 
to the OMs can be considered, which 
affect o"oM) e and o-^^™^*^: 

• coating the OMs with wavelength- 
shifters would increase the photon- 
collection efficiency without affect- 
ing the noise rates. 

• using PMTs with a larger cathode 
area would increase the efficiency, 
even though the noise level would 
go up. 

• reducing the potassium levels in 
the pressure glass would reduce 
the dark noise spread. 

• Fig. 3 shows that there are ice- 
depths with larger transparency 
which are most suitable for SN de- 
tection. In fact, the transparency 
is expected to be extremely good 
at shallower depths [37,38]. How- 
ever, the muon flux is larger closer 



to the ice surface, producing more 
background Chercnkov light. 
All these improvements have been 
investigated and we estimate that 
they can reduce the value of cr^™^^ by 
20 - 40% each. 

Detailed investigations on after- 
pulsing in the range up to 1 msec 
are currently being pursued in order 
to optimize the choice of artiflcial 
dead-time. The effects of raising the 
PMT thresholds remain to be inves- 
tigated. Finally, it may be noted that 
a sliding time-window, fltted to the 
supernova time-proflle, rather than 
applied at arbitrary times, would in- 
crease the galactic coverage of the 
detector studied in this paper by an 
estimated 11%. One could even op- 
timize the size of the window to flt 
the peaked part of the profile so that 
the ratio e/o"oM is maximized. An 
additional 10% could be gained in 
this way, even though it is clear that 
such optimization is partly model- 
dependent (for the model used here, 
a 4 sec time window would be opti- 
mal). 

An uncertainty in the strength of 
the signal is the Pg temperature. In 
this paper, we use the same assump- 
tion as in [11] {Tiy^ = 4 MeV), but 
there is no precise prediction for this 
value and it may be lower (see e.g. 
[39] ) . Furthermore, corrections to in- 
verse beta-decay cross-sections due 
to recoil and weak magnetism [40] 
would lower the resulting positron 
energy somewhat. On the other 
hand, this can be compensated in 
some measure by the performance 
of the present detector (680 OMs), 
which is expected to be higher than 
the configuration studied in this pa- 
per. With the planned transition to a 
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detector of kilometer-cubed size, us- 
ing thousands of OMs, full coverage 
of the Galaxy will be attained. 

One of the greatest assets of 
Amanda with its low background 
rate is that it has the potential to 
see a supernova independently from 
other detectors. In the event of a 
discovery, the time profile and inten- 
sity of the burst could be studied. 
These aspects will be investigated 
in the future. Should an event occur 
in our Galaxy, the information that 
could presently be provided to the 
rest of the supernova community, is 
an estimate of AyU and a time stamp 
associated with the event. 

9 Conclusions 

Amanda primarily aims at detect- 
ing high-energy neutrinos. However, 
by monitoring bursts of low-energy 
neutrinos, it is operating as a gravi- 
tational collapse detector as well. In 
the course of this analysis, a method 
has been developed and applied to 
filter the data and remove unstable 
periods and OMs. from disturbances. 
Furthermore, the x^-bascd analysis 
method used here takes into account 
the varying statistical properties 
of different OMs. The sets of data 
taken in 1997 and 1998, correspond- 
ing to 215 days of live-time after 
cleaning, have been fully searched 
for supernova bursts occurring in 
our galaxy. No significant candidates 
were found. A conservative model of 
the supernovae-distribution in space 
and time has been used to assess the 
performance of the detector, result- 
ing in an estimated coverage of the 
galaxy of 70% with 90% efficiency. 



with one expected background fake 
per year. This performance is satis- 
fying, the more so since the primary 
design of Amanda is for high-energy 
neutrino detection. A real-time su- 
pernova detection algorithm, has 
been developed and tested on the 296 
days of live time available. The al- 
gorithm turned out to be sufficiently 
robust and resulted in five fakes with 
the set of cuts chosen. We find that 
800 OMs of the same type as those 
used in the first four strings of the 
array would suffice to cover the en- 
tire galaxy. It is also shown that this 
number can be reduced significantly, 
either by reducing the spread of the 
dark noise, or by improving the col- 
lection sensitivity of the OMs. 
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